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Abstract 

We consider solutions of two-dimensional m x m systems hyperbolic conservation laws that are con- 
stant in time and along rays starting at the origin. The solutions are assumed to be small L°° perturba- 
tions of a constant state and entropy admissible, and the system is assumed to be non-strictly hyperbolic 
with eigenvalues of constant multiplicity. We show that such a solution, initially assumed bounded, must 
be a special function of bounded variation, and we determine the possible configuration of waves. As a 
corollary, we extend some regularity and uniqueness results for some one-dimensional Riemann problems. 

1 Introduction 

1.1 Outline and motivation 

Consider a Id system of conservation laws 



Ut + /(£/). = 



augmented with an entropy inequality 



v{U)t+^{U),<0 

for an entropy pair (77, -0) with convex 77. Consider Riemann problem initial data 



U{x,t = 0) 
The solution is expected to be self-similar. 



U-, x<0 
U+, x>0 



X 
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U{x,t) = U{0 , ^ 
which reduces the system to 

[f{u)-^u]^^u , mu) ~ omh < v{u) 

which, for smooth U, is equivalent to the non-divergence form 



In prior work [14| we assume: 

1 . U has values in a closed ball Ve C V oi small radius e around U € R"* . 

2. The system is strictly hyperbolic: for every U G Pe, fuiU) has m real and simple eigenvalues A"(J7) with 



*This material is based upon work partially supported by the National Science Foundation under Grant No. NSF DMS- 
1054115 and by a Sloan Foundation Research Fellowship. 
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associated right eigenvectors r°'{U). 

3. Each eigenvalue is either linearly degenerate^ 

VC/eT'e : A^(C/)r"([/) =0 , (2) 

or genuinely nonlinear. 

yu c.v.-.xtjiuy'^iu) . 

Under these conditions we show, among other resuhs, that 

1. U€ SBV C BV (as a function of 0, and 

2. for ^ « A" where A" is Unearly degenerate, U is constant on each side of a contact discontinuity. 

A key step of the proof was to show the following is impossible: for linearly degenerate A" and / a closed 
interval of positive length, 

VCe/:A"(;7(0)=e . (3) 
This is straightforward if we already know U e BV: then U is differcntiable in at least on(0 ^ G /, so 

Q^[fu{U{i))-mUdO 
^[fu{U{i))-r{Um]Um 



UdO II r'-iUiO) 



But taking of ([3]) yields a direct contradiction: 

>^uiu{mm = 1 

(There is no contradiction if / is a point (so that 9^ cannot be taken), which corresponds to the familiar 
case of a contact discontinuity.) 

The obstacle was to prove — not assume — that U G BV, merely starting from U G L°° which docs not 
imply differentiability. A key ingredient is a result of Saks [26|: for an arbitrary scalar g : / — > M there is a 
sequence (^„) C / converging to so that U\d, U restricted to I? = {^o} U S, is differcntiable in Un- 
important systems like the 2d steady non-isentropic Euler equations with supersonic background state 
U (where one space direction acts as time t and the other Appendix B) have a linearly degenerate 

eigenvalue of multiplicity p > 2 (physically corresponding to contact discontinuities that have, in addition to 
tangential velocity jumps, entropy jumps as well). Including passively transported y, z velocities yields even 
higher multiplicities p. While many results for strictly hyperbolic systems generalize to p > 2 in relatively 
straightforward ways, this is not the case here. Reason: Saks result is not true for g : I ^ M.P with p > 2. 



In this case, 12 1 shows C"'^ continuity of U\jy in with optimal exponent f3 = - < 1. 



1.2 Related work 

There has been some related work on determining the regularity of weak solutions. Though BV is the 
standard setting for one-dimensional systems of conservation laws, a result due to Ranch [2J] shows that 
BV is not appropriate for unsteady multi-dimensional problems. An example of a one-dimensional SBV d 
regularity result is [l| , in which Ambrosio and De Lellis showed that (not self-similar) L°° entropy solutions to 
scalar one-dimensional conservation laws are special functions of bounded variation (as a function of x, for all 
but countably many values of t). In Q, Dafermos showed that BV self-similar solutions to one-dimensional 
systems with genuinely nonlinear fields are SBV (without reference to entropy). In Bianchini and 
Caravenna showed that BV entropy solutions to one-dimensional systems with genuinely nonlinear fields are 
SBV for all but countably many values of t. Though we require self-similarity and a smallness assumption, 
the fact that we can treat systems with degenerate fields and only assume L°° is interesting. 

^in fact almost everywhere 

^ The Lcbcsgue decomposition expresses a BV function as the sum of an absolutely continuous function, a singular (dif- 
ferentiable almost everywhere with zero derivative) continuous function (referred to as the Cantor part), and a jump (saltus) 
function. A special function of bounded variation (SBV) has vanishing Cantor part. 
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1.3 Overview 



In this paper we extend the results in IJ] to non-strictly hyperbohc systems with constant multipUcity. In 



Appendix [B] we verify that the required assumptions are satisfied for the full Euler equations. Though the 
regularity of steady and self-similar solutions of full Euler was proved by the second author in 25| , that study 
was not limited to small perturbations. Therefore, there was less that could be said regarding the structure 
of solutions, since strong shocks and subsonic flows were allowed — there was no notion of "sectors" (defined 
in Section I12p in which we know all waves of a given family must occur. That analysis used the algebraic 
properties of the full Euler equations instead of the implicit function theorem, and in this paper we treat 
general systems as well as verify that the full Euler equations fit into this perturbative framework in order 
to get more structural information than was available in [25[ . 



2 Physical Systems and Entropy Solutions 

A system of two-dimensional conservation laws is a system of nonlinear partial differential equations of the 
form 

Ut + r{U)^ + P{U)y = Q. (4) 

The unknown U{t,x,y) = ([/^(t, x, y), [/^(i, x, y), f/™(i, x, y)) is a function from R+ x R'^ to P C M™, 
the individual components 11°', a = 1, ...,m are called the conserved quantities, the set V is called the phase 
space of physically allowed values, and the smooth functions /^,/^ : V — > M™ are called the horizontal 
and vertical flux functions, respectively. (Like U, f^ and arc each column vectors with components 
/-",p«,a- 1,...,™.) 

Smooth functions ri,Tp^ ,Tpy : V ^ M. are an entropy- entropy flux pair for the system if, for all U £ V, 

r/i/jj is positive definite, (5) 
=Vufu, and (6) 

ru^mfl- (7) 

Definition 1. A physical system is a choice of conserved quantities, phase space, flux functions, and entropy- 
entropy flux pair as described above. 

Definition 2. U £ L°°(M+ x R^;^') is a weak solution to @ if for any test function $ = ($i $2 ^ 

^tU + ^.J^{U) + ^yfy{U)d{t,x,y)^Q. (8) 

Unfortunately, weak solutions are not unique for given initial data. Therefore, we need an admissibility 
criterion for weak solutions. In this case, the appropriate requirement is that 

etTl{U) + Q..r{U) + Qyr{U) d{t,x,y)<0, (9) 

2 

for any non-negative test function 9 G Cj?°(R+ x M^;M). This is what is usually referred to as the integral 
form of the differential inequality 

7?(f/)t + ^^(f/), + 7A^([/), <0. (10) 

Definition 3. [/ e L°°(M+ x R^; R"') is an entropy solution to a physical system (|4]) if it is a weak solution 
that also satisfies © for all non-negative test functions Q. 

3 Steady and Self-Similar Solutions 

We are interested in entropy solutions that are steady in time and constant on rays emanating from the 
origin. Using the same arguments as in [l^ . we have the following. 
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Definition 4. A steady and self-similar entropy solution U £ L°° to a physical system Q satisfies, in the 
sense of distributions, 

{fy{u)~triu))^ + f-{u) = o, 

i^HU) -^4''^{U))^+^P^iU) <0, x>0 . (11) 

(Once either the right or left half plane is chosen, U is only a function of ^ = y/x.) 

Remark 1. We will later justify why we can ignore the case of a; = without loss of generality. 
For the remainder of this chapter, we assume the following. 

Assumption 1. The system of conservation laws under consideration, is a physical system, and all 
solutions considered are steady and self-similar entropy solutions. 



4 Smallness and Intuition 

Many of our results require the implicit function theorem, and therefore a smallness assumption. To that end 
we assume that our phase space "P is a small neighborhood of a constant background state U , and rename 

it Ve- 

Assumption 2. The phase space of allowed values for the conserved quantities is of the form 



\U- 



U\<e}, (12) 



Ve := |C/ e ] 

for some small e > 0. Thus the solutions we consider satisfy 

\\U{-)^U\\l'^ <e. (13) 
We will reduce e as necessary throughout, but only finitely many times. 

Recall that we are not assuming any regularity of the entropy solution [/(•) — only that it is bounded. 
Therefore we cannot expect it a priori to be differentiable anywhere. We are not assuming it is of bounded 
variation either, and so we cannot even analyze its derivative in the sense of measures (or even talk about 
left or right limits). However, if we look at the differentiated form of the equations anyway, we obtain from 
the first Une of ([TT]) 

{f!j-m)Ui^O. (14) 

Therefore, if on some interval ^ is not a generalized eigenvalue of the matrix pair (^fij{U{^)),flj(U{^))^, 

= and therefore U is constant. If instead there is some interval of ^ on which ^ is an eigenvalue, then 

must be the associated eigenvector, which is precisely the case of a simple wave for the steady problem. 
Either way, the smallness assumption on the phase space suggests that when ^ is not close to a generalized 
eigenvalue of the matrix pair (/(/(£/),/[/(£/)), U must be constant. 

Shocks and contact discontinuities are also possible, but the smallness assumption and standard facts 
about conservation laws suggests that these waves also occur near the eigenvalues evaluated at the background 
state U, since the phase space is a small neighborhood of U. 

All together, the differentiated form suggests that all interesting behavior must occur near the eigenvalues 
of the background state. So even if we were allowed to use the differentiated form it would be important to 
analyze the characteristic behavior of our system, which is the next step. 



5 Pointwise Information 

It is cumbersome to work with the integral form of a conservation law, so we instead investigate what 
pointwise information we can derive from it. Recall that we are assuming that U £ L°°. Rearranging the 
first line of pT|) yields 

(/^(c/)-er(c/))^ = -r(t/). (15) 
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Therefore, the quantity 



(16) 



has a distributional derivative that is L°° (since is smooth on V), and is therefore almost everywhere 
equal to a Lipschitz continuous function. Therefore, the fundamental theorem of calculus asserts that 



r{U{Tj)) dv, fora.e. 6,6 
Similarly, the second line of PT|) shows that the distributional derivative of 



(17) 



(18) 



is a non-positive distribution, and therefore a non-positive measure. Therefore, there is a version of psp 
that is a non-increasing function of bounded variation, and rearranging this statement we obtain 



< ilj'^{U{r/)) dr], for a; > and a.e. ^ < ^ 

hi 

Similarly, we can obtain 



(19) 



rS.2 

> ip''{U{r])) drj, for .t < and a.e. ^ < 6. 



(20) 



We now fix a version of U so that these equations and inequalities hold for all ^ < ■^2, as discussed in 
Section 7 of 141, so that 



(fy{u)^cnu)) 

{r{u)-ir{u)) 



?2 



«2 



«2 



r(t/(^)) din, 



?2 



< - / V{U{ri)) drj, x>0, 
> - r r{U{ri)) dTj, x<0, 



(21) 



for all 6 < 6. 



6 Hyperbolicity 

Consider the homogeneous polynomial 

Fix : y) det (f x fu{U)) = det (x/^(C7) - yf^{U)), 



(22) 



where (x : y) are homogeneous coordinates on RP^, x = {x,y), and / = (./^,/^)- We call the background 
state ?7 hyperbolic if P(a; : y) has exactly m roots in KP^, counting multiplicity. If i? is any rotation matrix, 
it is easy to see that Rx x Rfu ~ xx fjj. Therefore, if / is replaced by Rf, the roots of P are rotated by the 
same amount. P has degree at most m, so there are at most m distinct roots and therefore we can rotate / 
to ensure that (0 : 1) is not a root of P. Assume without loss of generality that this has been done. 
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The fact that (0 : 1) is not a root of P immediately imphes that 



det/^^(C/)#0. (23) 

Now consider the polynomial 

p(0:=det(/^(I7)-C/^(I7)). (24) 

From the above discussion, p has m real roots, counting multiplicity, since each of the m roots of P lead to 
a finite (since none lie on the y-axis) value of ^ that is a root of p. A generalized eigenvalue X{U) and an 
associated generalized eigenvector r{U) of the matrix pair {fij{U),fif{U)) satisfy 

{fy{U) - X{U)f^{U))r{U) = 0. (25) 

(Note that the term generalized eigenvectors in this context refers to the elements of the kernel of the linear 
matrix pencil {—\f^ + /^), not elements of the kernel of (— A/^ + f^,)'^ for fc > 1 in the context of defective 
geometric multiplicity.) 

For the remainder, we shall assume that the steady problem is hyperbolic in the following sense. 

Definition 5. The steady problem associated to the system Q is called hyperbolic on the phase space 
Vc if the generalized eigenvalues of the matrix pair {fij{U), fjjiU)^ are real, semisimple, and of constant 
multiplicity on (thus the sum of the multiplicities equals m). It is called strictly hyperbolic if there are 
m distinct generalized eigenvalues, and non-strictly hyperbolic with constant multiplicity otherwise. 

Assumption 3. The steady problem associated to ^ is non-strictly hyperbolic with constant multiplicity 
on the phase space (which implies flj{U) is non- degenerate for all U G Ve)- 



As discussed before, the strictly hyperbolic case was treated in |14| . and the purpose of this paper is to 
treat repeated eigenvalues. 

7 Change of Dependent Variables 



In light of Assumption [21 we can perform a change of dependent variables as in |14| . Since /^(J7) is non- 
degenerate, 

U^V:=f^{U) (26) 
is a diffeomorphism on V^, after reducing e if necessary. We then define 

V:^f^{U), f{V):^fy{U{V)), (27) 

eiV):=r{UiV)), qiV):=^py{UiV)). (28) 

Then, 



Also, we have 



and 



fv ^flUv. (29) 



ev = ^uUv = VufuUv = Vu, (30) 



qy = tp'fjUv = VufijUv = evfv- (31) 

Therefore, e and q form an "entropy-entropy flux pair" for the flux /. The term is applied loosely here 
because e is not necessarily convex. Properties of the entropy are only needed in two instances, and further 
properties of e will be discussed when they are needed. 

Abusing notation, we shall continue to refer to our phase space as Ve, but it will now refer to a small 
ball around V of permissible values. 
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The differential equations and entropy inequalities for V are then 

f {fiv)~^v)^ + v = o, 

{qiV)-^eiV))^ + eiV)<0, x>0 
{q{V)~^e{V))^ + e{V)>0, x<0 



(32) 



The pointwise form is 



(fiV) 


SI 


iiiv) - 




('Z(^) - 


-Ce{V)) 



«2 



< - 



> 



V{ri) drj, 

?2 



e(F(77)) drj, x>0, 



(33) 



for all ^1 < ^2- 

8 Eigenvalues and Eigenvectors 

Since fv = f!jUv = W^)'', 

det(/^ ~ A/^) - ^ (dct(/f,-A/^))(dct(/^)-i) =0 ^ dct(/y - A/) = 0. 



(34) 



Therefore, the generalized eigenvalues of the matrix pair {f^, ffi) are precisely the eigenvalues of the matrix 



To that end, define 



to be the distinct values of A solving 



X\V) < X'^iV) < ■ ■ ■< A"(y) 



(35) 
(36) 

(37) 

the eigenvectors 



det(/v(F)-A/)) -0. 
If s is a generalized eigenvector with eigenvalue A, then 

= (/^ - A/^)s - - Xf5){f&)-'f^s ^ {fv - A/)/^s, 

which means that f^s is an eigenvector of fv- Since the generalized eigenvectors span 
of fv do as well. Define 

R"{V) kcr {fv{V) - X''{V)l), p„(F) dim i?"(F), (38) 

so that, under the hypcrbolicity assumption 

p^(V)=:pc.onP„ (39) 

and 

n 

R" = R"{V) for aU V eVe- (40) 

Q = l 

In the strictly hyperbolic setting, it is relatively easy to prove that if the matrix fv is a smooth function 
of V, then so arc the eigenvalues and eigenvectors (see for example [Hj). However, the situation is more 
delicate in the case of repeated eigenvalues — there are examples in which the eigenvalues and eigenvectors 
arc not as smooth as the matrix. Fortunately, if the eigenvalues are semisimple and of constant multiplicity, 
they and the eigenvectors can be shown to be as smooth as the matrix, though the individual eigenvectors 
are only guaranteed to be locally defined as smooth functions (the eigenspaces are smooth when given a 
suitable topology). Since we are only considering small perturbations, we can simply reduce e if necessary 
and have our right and left eigenvectors defined on all of Ve- 



7 



The smoothness of the fluxes, the hyperbohcity assumption, and the discussion in Section in the 
Appendix aUow us to conclude that, for each a = 1, .., n, 

A" : n ^ K (41) 

is smooth. In addition, we have for all V € Ve an orthonormal basis for 

=Span {r"^i(F),...,r"'P°(F)}. (42) 

Reducing e as necessary, we have that the right and left eigenvectors 

r°'''{V),r^'{V) -.Ve^W (43) 

are smooth, and satisfy the normalization 

\r°''\V)\=l, (44) 
r'\Vy'^iV)=d^0S,, Va,/3 = l,..,n,z = l,..,p„,j = l,..,p^. (45) 

If for some a, = 1, then we omit the second index of the eigenvector and simply denote it by r". 

9 Genuine Nonlinearity and Linear Degeneracy, Convexity of e 

The results in this paper are systems with eigenvalues that are either linearly degenerate on all of Ve or 
genuinely nonlinear on all of V,,- 

Definition 6. An eigenvalue A" is linearly degenerate if 

XviV)r°'^'{V)=0 onVe, i = l,..,Pc.. (46) 

As it turns out, if Pa > 2, A" must be linearly degenerate, and there is a nice geometrical structure 
created by the eigenspaces. This result is originally due to Boillat, and this theorem and its proof can be 
found in [27|. 

Theorem 1 (Boillat, as in (27j'). Suppose the hyperbolicity assumption (Assumption\^ is satisfied. If an 
eigenvalue A" has multiplicity Pq > 2, then it must be linearly degenerate. In addition, the affine subspaces 
V + R"{V) are the tangent spaces to a family of sub-manifolds of dimension p^. Each integral submanifold 
can be parameterized by s W"'{V~' , s), with s in M^°, so that 

W°'{V~,0) = V~ (47) 
Span [w^,{V-,s)Y'' ^ = R°'{V-). (48) 

They form a foliation of Ve called the characteristic foliation associated with A". 

We now recall the definition of a genuinely nonlinear eigenvalue, which by the previous theorem must 
have multiplicity 1. 

Definition 7. An eigenvalue A" is genuinely nonlinear if 

XviV)r°'{V) ^0 onVe- (49) 
We can orient r°'{V) so that, without loss of generality, 

Ay(F)r"(V^) > on-pe- (50) 

We make the following assumption. 

Assumption 4. Each simple eigenvalue is either genuinely nonlinear or linearly degenerate. {Recall that 
an eigenvalue with multiplicity greater than one must be linearly degenerate, so no assumption is necessary 
in that case.) 
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We will, as in need to construct an averaged matrix to continue the analysis. To ensure the existence 
of an appropriate averaged matrix, we need the "entropy" e to be convex (or concave so that — e is convex). 
Since e is the original horizontal entropy flux composed with the inverse of the horizontal flux, and it is the 
original entropy 77 that is convex, wc must carefully inspect whether e is convex or not. 

Lemma 1. If f^{U) has only positive (negative) eigenvalues, then e is strictly convex (concave). 

Proof. We shall use Proposition 6.1 from It states that if H is symmetric positive definite, and K is 
symmetric, then HK is diagonalizable with real eigenvalues. Moreover, the number of positive (negative) 
eigenvalues of K equals the number of positive (negative) eigenvalues of HK. First, from (PO)) . we have that 

evv=riuuUv, (51) 

which we rewrite as 

ifSr' = (vuu)~'ew (52) 

By assumption, ('i]uu)~^ is symmetric positive definite, and eyv is symmetric. Then, applying the proposi- 
tion, if aU the eigenvalues of (f^)~^ are positive (negative), then eyv is positive (negative) definite, since a 
symmetric matrix is positive (negative) definite if and only if its eigenvalues are all positive (negative) . □ 

We finally note that changing variables to V docs not affect linear degeneracy or genuine nonlinearity. As 
we showed before, if s is a generalized eigenvector of the matrix pair (/^, /^), then r = f^s is an eigenvector 
of /y, and so 

\yr = \v(f^s) = \uUvVus^\us. (53) 

This docs not matter much when discussing general systems, but when checking for linear degeneracy or 
genuine nonlinearity in a specific system it is more natural to check in terms of the original conserved 
quantities U or perhaps another choice of state variables. 

10 Averaged Matrix 

To proceed with the analysis, we need to construct an averaged matrix A that is 

• smooth and diagonalizable (with real eigenvalues), 

• satisfies A(V- ,V+)(V+ - V') = f(V+) - f(V-), and 

• satisfies A(V,V) = fviV). 

A commonly used choice that suffices in [3l is to use 

A(V-,V+):^ f fv{V- +s(V+ -V'))ds. (54) 
Jo 

Clearly this satisfies most of the requirements — but in fact this definition only guarantees a diagonalizable 
A in the strictly hyperbolic case. This is because the set of matrices with all simple real eigenvalues is open, 
and wc are only interested in G Pc, and so smoothness of the flux guarantees that this averaged matrix 
is a small perturbation of fv(V). However, the set of matrices with repeated real eigenvalues is not open, 
and so in the repeated eigenvalue case this A is not guaranteed to be diagonalizable. 

An averaged matrix of this type is often used in numerical computations, and for some specific systems 
there is a Roe averaged matrix available. It has the special property that it can be defined as 

A(V-,V+) = fv(V), (55) 

where V is some appropriate averaging of the states V~ and In this case, diagonalizability is guaranteed 
from diagonalizability of fv- (It also has the useful property that expressions for the eigenvalues and 
eigenvectors of fv are available, which makes it especially well suited for numerical computations.) This 
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is often accomplished by doing a line integral in phase space between the two states, but choosing a more 
sophisticated path than simply the line segment between the two states which allows one to analytically 
evaluate the integral. However this is not available for general systems since it is reliant on the algebraic 
structure of the flux functions, and so we need another approach. 

A theorem due to Harten and Lax that appeared in [18| shows that physical systems always possess an 
averaged matrix with these properties. Therefore, we need either ew to be positive definite, or negative 
definite so that (— e) can function as a convex entropy. In light of Lemma [1] we will need to assume that 
the eigenvalues of f^{U) are all the same sign (from which it follows that they will all have the same sign 
on all of Ve)- 



Theorem 2 (Harten, Lax as in [18[). Suppose there exists an entropy /entropy- flux pair {e,q) for the flux 
function f with eyy positive definite. Then we can define an averaged matrix A{V^), such that it is smooth 
in , A{V,V) ~ fv{V), o,nd it is diagonalizable with real eigenvalues for all G Ve- Most importantly, 

fiV+) - fiV~) = AiV-,V+){V+ - V-). (56) 

A key idea in the proof is to use a real symmetric square root of the matrix evv i which is only possible 
if evv is symmetric positive definite. Since we are assuming that at least one eigenvalue has multiplicity 
greater than one, we need to make the following assumption in order to have an averaged matrix. 

Assumption 5. Assume that the eigenvalues of f^(U) are all positive or all negative. 

Wc denote the eigenvalues of A{V^) as A"''. If for some a, pa = 1, then from the discussion in Section 1X1 
in the Appendix it follows that A" and f" are smooth functions of V^. If instead Pq, > 1, we will not have 
in general that A"'^ = ... = A"'P°. Since the multiplicity of these eigenvalues is not necessarily constant on 
Ve X Ve, we cannot conclude that these eigenvalues and their associated eigenvectors are smooth functions 

of y±. 

However, as discussed in Section [A] in the Appendix, the eigenvalues are continuous functions of . We 
define, for a — 1, ..,n, the a-group to be 

{A"'\A"'2,..,A"^f°} , (57) 

which can be continuously labeled so that 

X°''\V,V) = X"''^{V,V) = ... = X"''P''{V,V) ^ X"{V) (58) 

for all V eVe- 

Many examples exist that demonstrate lack of continuity of eigenvectors when the multiplicity of an 
eigenvalue changes, even for symmetric matrices. Moreover, the projection operators onto the eigenspaces 
also can display this behavior — it is worse than just not being able to smoothly pick bases for these 



eigenspaces. However, if we instead consider, as in [22|, the total projection for the a— group, then this will 
be as smooth as A. We have the formula 

PT{V^)^^JjyzI-A{V^))-^dz, (59) 

where Pp is the sum of the projections onto the eigenspaces of all eigenvalues inside some counterclockwise 
contour F. By continuity, all eigenvalues in the a-group remain close to X°'{V) for all in Ve, and so we 
have for a = 1 , . . , n that 

P"[y^)^^- I _ {zl ~ A[V^))-^dz (60) 



27ri 



\z~\^(V)\=5 



is the total projection for the a-group, where 5 is small enough so these curves remain distinct for different 
choices of a, and large enough so as to include the entire a-group for all 6 Ve- Clearly such a 5 > can 
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be found for e sufficiently small. We will make use of the following properties of the projections. 

P°'pP=SapP", (61) 
= iP" = P'^AP", (62) 



(63) 



Note that this implies that P"M'" is an invariant subspace for A. 

If for some a, Pct = 1, then A" and its associated right eigenvector are smooth functions of V^. In 
this case we also normalize so that 

f°'{V,V) =r°'{V), (64) 
|f"(l/^)| = 1, for aU a with = 1. (65) 

Under the assumption [SJ we can rewrite the pointwise version as (after slight rearrangement) 



{q{V)~Ce{V)) 
(qiV)-CeiV)) 



i2 



< - / e{V{rj)) drj, 
> - / e(F(77)) drj, 



X > 0, 
X <0, 



(66) 



for all ^1 < ^2 



11 Left and Right Sequences 

We now continue as in [l3|. Since V{-) is only assumed to be bounded, it is not necessarily piecewise smooth, 
and and at each point it is not necessarily either continuous, approximately continuous, or discontinuous 
with an approximate jump discontinuity (these are the hypotheses typically used in deriving the familiar 
Rankine-Hugoniot jump conditions). Arbitrary bounded functions do not even possess well defined limits 
from the left and right, and so we have to be more careful in this setting. 

We wish to use the pointwise form (|66p . so to that end we define pairs of sequences (^^), {£,^) such that 
C < ik for all k and |^ 

Since Ve is compact, there exists a subsequence ^^^^ such that V{£,ffi.^) V'^. Similarly there exists 



a subsequence such that l^(^iQ(fc))) 

sequences such that 



^ V . By construction, taking 
F(4±) V^, 47 < for all k. 



(67) 



Understanding that this depends on the pair of sequences chosen, we define 

[g{V)] :^ g{V+) - g{V-) 

for any function g V. Define 

J(.g(V^);0 :=sup|[<7(V^)]|, 
where the supremum is taken over all such sec^uences. Obviously J{g{V); ^) = 

Applying ([55)1 with ^ ~ ^ and ^2 ~ (this is why we insisted 4, < for all k) and taking the limit 
fc — > 00 yields 

[q{V)] - ^[eiV)] < 0, x > 0, (70) 

[q{V)]-^[eiV)]>0, x<0, 



(69) 

g oV is continuous at 
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with the first hne being equivalent to 

inv)] - av] = 0, 



(71) 



which is the famihar Rankine-Hugoniot condition. Therefore, even in the absence of well defined left and 
right limits, we can make some sense of limits using these pairs of sequences, and the usual conditions apply 
to these sequence-dependent V^. 

The first line of (|66p implies that, for a given pair of sequences 

[V] =0 or C = A"'* and [V] £ ker(i(F±) - A"^V). (72) 



12 Sectors 

We now start to make rigorous statements similar to the intuition developed in Sectional and present two 
theorems from [iJl relying on the existence of A. 



Theorem 3 (EUing, Roberts [IJ]). Suppose V is continuous on an interval I =]^i,^2[ o,nd that ^ is not an 
eigenvalue of fv(V{^)^ for any ^ S /. Then V is constant on I. 

The continuity assumption is stronger than we need, so we proceed without assuming continuity but 
instead assuming ^ is bounded uniformly away from all eigenvalues of fv{V{^)). 



Theorem 4 (Elling, Roberts [lj|). Consider an interval I = ]^i,^2[- There is a Ss ^ '5s(e) > 0, with 

6s iO as eiO, (73) 

so that 

Va = l,..n,VC e / : |A"(F) -^1 >(5, (74) 

implies V is constant on I. 

Recall that under Assumption |3] 

P(0:l) = det/^(I7)7^0. (75) 

Find some ^ such that 

P(l:C)=det(/^(C7)-C/f5(C7))7^0. (76) 

Rotate coordinates so that (1 : ^) becomes (0 : 1) and (0 : 1) becomes (—1,^). In these new coordi- 
nates. Assumption |3] is satisfied, and so all results apply. Since P is invariant under rotation, in these new 
coordinates 

P(-l,O7^0, (77) 

which means that — ^ ^ A"(y) for all a = 1, .., n. Reduce e so that | — ^ — A"(V^)| > Ss for all a. This will 
still be true on some small interval around — ^, and so the previous theorem applies. Rotate back to the 
original coordinates to deduce that V must be constant in thin sectors containing the positive and negative 
y-axes. Therefore, we lost no generality in only considering test functions supported away from a; = and 
doing all the analysis in terms of 

Now, construct n intervals of the form 

/":= ]A"(y)-5„A"(y) + <5,[. (78) 

y 

Considering both x > and a; < 0, we now have 2n thin sectors centered at — = X"{V) for some a. V 

X 

is constant outside these sectors, and we have 

lim 1/(0, - lim ViO- (79) 

x>Q Jc<0 

All interesting behavior occurs within these sectors, so we will consider the behavior of V in these sectors 
individually. 
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13 Linearly Degenerate Sectors 
13.1 Satisfying the Jump Conditions 

We now analyze the possible behavior of V in I", where A" is linearly degenerate. 

The first thing we recall is that if V is discontinuous at ^, then we can find a pair of subsequences whose 
limits yield 

IfiV)] - m = 0, (80) 

with [V] ^ 0. We need to show that ([50]) is satisfied if and only if both lie on the same leaf of the 
foliation discussed earlier. 

First, we notice that the eigenvalue A" is constant on each leaf of the foliation. By direct calculation, we 
see that 

A(a"(M^"(F-,s))) = X^{W"{V~,s))w:.{V-,s) ^ (81) 

for all i e {1, ...,pc(} by linear degeneracy (since £ i?"). Therefore, 

\"'(W°'{V^,s)) is constant. (82) 

We claim that for all s, if V+ = W"{V-,s) and ^ ^ \°'{V-) = \°'{V+) then §Q is satisfied. Make 
these choices for and ^ and consider 

F{s) :=/(Ty"(y-,s)) -/(y-)-A"(M/"(V^-,s))(M/"(y-,s)-y-). (83) 

Notice F(0) = 0, and (using ((5T|) ) 

Fs^ is) = fv [V- , s))w:. (W" {V- , s)) - A" (ly" [V- , s))w:^ (y- , s)) ee 0. (84) 

Now we check for entropy admissibility. Define 

E{s) ■~q{W°'{V^,s)) -q(V~)- X°'{W°'{V-,s))(^e{W°'{V-,s)) -6(1/")). (85) 

Then £'(0) and 

Es^ (s) = qv{W"{V-,s))W^,{V-,s) ~ X''{W''{V-,s))ev{W''{V-,s))W^,{V-,s) (86) 

= ev{W''{V-,s))(^fv{W"{V-,s))~X"{W"{V-,s))l)w:.{V-,s)=0, (87) 

(using the property of entropy-entropy flux pairs). Therefore, this choice of and ^ satisfles for 
either x > or a; < 0. This means that any two states on a leaf of the foliation can be the left and right 
sides of a contact discontinuity located at ^, if ^ is the (constant) value of A" on that leaf. To that end, we 
will sometimes refer to W°'{V~ .^)sls the contact manifold through V~ . 

We now use a theorem from [l^ to prove that these are the only choices that satisfy ([50)) for e sufficiently 
small. 

Theorem 5 (Freistiihler [l^). For any (y~ in a sufficiently small neighborhood of {V , V , X" (V)) , 
if (V~ satisfy the Rankine- Hug oniot jump condition then V~ and V'^ must lie on the same leaf of 

the characteristic foliation, and ^ — A"(V^^) — A"(F+). 

We then have the following lemma. 

Lemma 2. For all ^ £ ^ i-> A"(F(^)) is continuous. 

In addition, ifi^ A"(V"(C)) on some open interval in I" , then V is constant on this interval. 

Proof. Suppose V is discontinuous at ^ G Then we can find a pair of subsequences so that [V] ^ 
and ([5n|) is satisfied. However, Theorem [S] applies and so lie on the same leaf of the foliation, and so 
[A"(F)] ~ 0. This holds for any pair of subsequences, and so A(T^(^)) is continuous at all S 

Similarly, if C 7^ then by Theorem [Sj (|80l) cannot be satisfied for any pair of sequences unless 

[V] = 0. Therefore, V itself must be continuous on such an open interval, and Theorem [3] shows it is 
constant. □ 
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13.2 Intermediate State 

Lemma 3. Let W''{V^ , s) be as above. For every V~ ,V G V^, there exists a unique s ~ s{V~ ,V) G 
-85(0) C such that 

P°'{W°'{V^,s), V) {V - s)) = 0, (88) 

(for some S > 0). 

Proof. Recall that P" has rank pa , so wc can view the map F defined as 

F{V-,s,V) -.^ P"{W°'(V-',s),V){V ~W"{V-,s)) (89) 

mapping ]RPa+2"> to W" . (More concretely, we may define F to simply be the Pa entries of F corresponding 
to the Pa linearly independent rows of P" , which do not change on Ve since linear independence is an open 
condition and e can be decreased.) Then, 

F,(y,0,7) = -P''(V,V)(w^^{V,0)\...\WXiV,0)). (90) 

However, recall that each W^i{V,0) is an eigenvector of fv{V), and therefore lies in the total eigenspace 
(which at {V, V) is just the eigenspace) that P" is projecting onto. Therefore, 

F, (F, 0, F) = - (w:^ (F, 0) I ... I {V,0)). (91 ) 

By construction of the contact manifold, this has rankpa, since the span of the columns is precisely R'^{V). 
By the implicit function theorem we obtain s{V~,V) close to the origin satisfying F{V~ , s,V) = for e 
sufficiently small. □ 

13.3 Regularity of /3 Components 

The following lemma establishes the regularity of the total projections along the other groups (3 a. 

Lemma 4. There exists a set E of full measure in I" such that for all S,o E, S, ^ I", P ^ 01, and 
i = 1, ...tPp we have 



{w'^{v{^,),s{o),v{o) (nc)-w^"(neo),s(e)) 



(92) 

(ie-eoi) + o(mo-^(eo)i-ie-Coi), 



provided that ^0 = '^"(^(^o)); o-nd defining s(^) := s(l^(^o)i ^(0) context of the previous lemma. 

Proof. From (|66p . we have 



(/(^(O) - /(^(Co))) - eo(^(e) - ^(Co)) = / V{i) - V{r,)dv. (93) 
We claim that the left hand side is equal to 



{f{ViO)-f{W'^{Vi^o),siO)))-^o{viO-W'^{Vi^o),siO)). (94) 

This is by design, since a contact at ^0 satisfies the Rankine-Hugoniot condition for the states V^(Co) and 
W°'{V{^o),s{0)- For ease of reading, abbreviate V{^o) = Vq, V{0 = V, and {V- , s{£,)) as W. Then 

{f{V)~f{V^))~UV^V^) (95) 

= {f{V) - f{W) + f{W) - fiVo)) - CaiV -W + W-Vo) (96) 

= (fiV) - fiW)) - Co{V -W) + {f{W) ~ fiVo)) - ^oiW - Vo) (97) 

^{f{V)^f{W))^^oiV-W) + 0, (98) 
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since = ^"(^0) by assumption. Then, using the averaged matrix we have 

(^A{W,V)~^o){V~W)^ f\~V{v)dii. (99) 



«0 



Left multiply both sides by P^{W, V) and add and subtract Vq inside the integral on the right to obtain 

{W,V) (^A{W,V) -^o'j{V -W) (100) 

= P^{W,V) I f^Vo-V{ii)dTj+{V-Vo)i(-^o)] (101) 



Lebesgue's differentiation theorem asserts that the integral on the right side is o(|^ — ^o|) for £,0 ^ E, a. 
set of full measure. Using the properties of the total projection we obtain 

(i(W^, V) - £o)P^{W, V){W - 1/) = o(|e - Col) + o{\v{0 - V{£,)\ ■ ic - Col) (102) 

However, since £q S and any eigenvalues in any /3-group are thus uniformly bounded away from 
(A(W, V) — Co) is uniformly non-degenerate on P^M™, and so 

(^A{W,V)-£,o^P'^{W,V){W~V) >SP'^{W,V){W-V), (103) 

for some S > 0, and the result follows. □ 

13.4 Regularity of V on Subsequences 

The following lemma applies the main result of [T^] to obtain a subsequence on which X^(C) is Holder 
continuous, although for our purposes continuity would be enough. 

Lemma 5. For almost every Co E E C I" , there exists a subsequence (Cn) Co o.'i^d < C < 00 such that 

m6o-^(Co)i <cic„-cor/". (104) 

That is, for almost all £0 £ E d I" there exists a set D :— {Co} U {(Cn)} such that V\d is Holder 
continuous with exponent 

Proof. As y(C) is a function from M to M™. the result follows directly from [l2|- □ 
Lemma 6. For (Cn) — > Co above, we have the following estimate for all (3 ^ a: 

p'3(i^"(y(Co),s(c„)),y(Cn))(^(Cn)-M/"(y(Co),s(c„))) =o(|c„-Col). (i05) 

Proof. The previous lemma proves that V\d is continuous, and so |l^(Cn) ~ ^(Co)| is o(l) and the result 
follows immediately from (IMI) . 

□ 

Lemma 7. In the setting of the previous lemmas, 

(y(C)-M^"(v^(Co),s(c„))) =o(|Cn-Co|). (106) 

Proof. We have 

(ncn)-w^"(neo),s(c„))) (107) 

= ^p'^(iy"(y(Co),s(c„)),y(Cn)) {viU - vi^"(^(Co),s(c„))) (108) 

= Y.P^{^"i^i^o),s{U),V{U){v{U-W"{V{£o),s{U)) (109) 

= o(|Cn"Co|), (110) 

due to the clever choice of s(Cn)- □ 
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13.5 Existence of at Most One Contact 



We now combine the results of the previous three sections to obtain the main result. 



Theorem 6. On a linearly degenerate sector, V is either constant, or constant on each side of a single 
contact discontinuity. 

Proof. By LemmaEl F := e I" \ £, = A"(y(C))} is closed and V is constant on 7"\F. 

Assume there are ^1,^2 ^ F and rj G I" \ F with ^1 < 77 < ^2- Then we can choose a maximal ,ri^[ 
containing jy but not meeting F. Necessarily 77^ G F, so 77+ = A"(y(77+)) and 77^ = A"(y(77^)). But V is 
constant on ]77~,77+[, so 77+ = 7/~, which is a contradiction. 

Hence F must be a closed interval. 

Assume F has positive length, and pick (,o G F with F D (^„) ^0 such that (|104p implies (|105p . Then, 



X"{V{^n))^X"{w''{Vi^0),s{^n)))+o(\vi^n)-W"{Vi^0),si^n))\) 



= eo + o(ic„-eoi). 

However, since ^„ G F, this implies 

- Co = 0(|C„ - Col), 

a contradiction. 

Thus F must be a point, or empty (which can be ruled out, but this is unnecessary). 



(Ill) 
(112) 
(113) 

(114) 
□ 



14 Genuinely Nonlinear Sectors and Global SBV Regularity 

Outside linearly degenerate sectors, the results in [l^ apply with the same proofs as in the strictly hyperbolic 



case, and so we have the following (see [IJ] for detailed discussions of shocks and simple waves) 



Theorem 7. Under Assumptions 1-5, and supposing that fij(U) has only positive eigenvalues, we have: 
(a) V must be constant outside of 2m thin sectors centered at the roots of (|24p . which we can group as m 
forward {x > 0) and m backward (x < 0) sectors. 

(h) Linearly degenerate sectors each contain at most one contact discontinuity. 

(c) Genuinely nonlinear forward sectors each contain at most one shock or simple wave. 

(d) Genuinely nonlinear backward sectors can each contain infinitely many shocks and simple waves, but 
there cannot be consecutive simple waves. 

(e) Each shock wave has a neighborhood on each side on which U is constant, and the size of the neighborhood 
is lower bounded proportionally to shock strength. 

(f) The width of the sectors and constant of proportionality depend only on the system and e. 

Remark 2. If fij{U) has only negative eigenvalues, then the "forward sectors" (named since they behave 
like forward-in-time Riemann problem solutions, in which at most one wave from each family appears) lie 
in X < while the "backward sectors" lie in a; > 0. 

Remark 3. If e is convex, we have shown that Assumptions 1-5 show that the solution in a; > must be the 
unique Lax solutior^ to the (one-spatial dimension y with x functioning as time) Riemann problem with data 
prescribed on the y-axis, and if e is concave, then this is true in a: < 0. Note that uniqueness is not expected 
in the half-plane corresponding to the backward-in-time solution, but our analysis shows the solution will 
still be of bounded variation. Since ([5^ are the same equations that a self-similar (function of C = x/t) 
solution to a one-dimensional system of conservation laws, our regularity results (and the uniqueness implied 
by our structural results and the uniqueness of Lax's Riemann solution ) apply to one-dimensional systems. 
We have therefore extended uniqueness from the class of self-similar solutions consisting of up to n waves 



^ Lax's solution is a forward-in-time self-similar solution for sufficiently close Riemann states consisting of up to n waves 
interspersed with constant states, and is unique among solutions having this structure. See jSJ, Section 4.6 for details. 
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with constant states in between to the class consisting of small self-similar perturbations of a constant 
state. This extends an earlier result of Heibig [lij . which considered forward-in-time Riemann solutions to 
onc-dimcnsional strictly hyperbolic genuinely nonlinear systems, and so our analysis can treat systems that 
can have repeated linearly degenerate eigenvalues. 

Proof. The statement regarding degenerate sectors follows from Theorem |6l and the remaining statements 
follow directly from [l^. □ 

In addition, the regularity results contained in [l^ also apply, now that we have at most one contact in 
degenerate sectors. 

Theorem 8. Suppose V satisfies the hypotheses of the previous theorem. Then V is of bounded variation. 
More specifically, V = Vs + Vl where Vs is a saltus (jump) function of hounded variation, and Vl is Lipschitz, 
and so V is a special function of hounded variation. The total variation of Vs and the Lipschitz constant of 
Vl are independent of V and only depend on the system and e. 

Proof. Again, now that linearly degenerate sectors have been treated by Theorem [6l and the analysis in [l3| 
works the same everywhere else, the claim follows. □ 



15 Comparison to the Strictly Hyperbolic Case 



Though much of the analysis is similar to the strictly hyperbolic case treated in |14| , there are some crucial 
differences. 

As discussed in Section [TOl the construction of the averaged matrix A{V^) is less straightforward in this 
non-strictly hyperbolic case and relics on convexity of e(-) = {{,f^)~^ {■)) ■ In [l3|, this assumption was 
not needed to construct the averaged matrix (the more straightforward definition (|54|) could be used), it was 
only necessary that 

eyyr"r" / (115) 

for each a corresponding to a genuinely nonlinear field. This quantity cannot vanish as long as is non- 
dcgcncratc, but if has eigenvalues with different signs then this quantity will have different signs for 
different fields. This was needed to identify which part of the Hugoniot locus corresponded to admissible 
shocks for genuinely nonlinear fields, and this ultimately determined whether the "forward" and "backward" 
sectors corresponded to a; > or a; < 0. Therefore, there are some strictly hyperbolic systems and background 
states that can be treated by Q for which neither x > or a; < (or any half-plane) is analogous to a 
forward-in-time Riemann problem solution to a one-dimensional system of conservation laws (recall that 
forward-in-time Riemann problem solutions consist of at most one wave from each family — which is the 
behavior described in Theorem [7] for forward sectors). However, since in the non-strictly hyperbolic case we 
need the eigenvalues oi f'fj{U) to all have the same sign in order to construct A, the systems and background 
states treated in this paper have the property that either x > or a; < will correspond to a forward-in-time 
Riemann problem solution. 

More importantly, the entire approach for linearly degenerate sectors is completely different in this non- 
strictly hyperbolic case. In the strictly hyperbolic case in the situation was similar to that in Theorem 

El 

X"{V{C}) = C on a closed interval F C and (116) 
P^{V{^o),V{0){V{0-V{^^))=0{\i^^^\) for ah a. (117) 

The approach was then to use a result due to Saks [2^: for any function F : M D A — > M, the set of points 
X at which 

AF{x + h)- F{x)\ , , 

liminf ' ^ ^ ^ ^ oo (118) 

h^o+ h 

is of measure zero. This implies differentiability at almost every S^q ^ F after restriction to a subsequence 
converging to of the real valued function 

e^p(y)(i/(0-neo)). (119) 
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Taking further and further subsequences for the to — 1 other (Lipschitz from (|117|) ) components allowed for 
the construction of a sequence (^„) — !• so that 

lim miLllM =: V' (120) 

exists and is finite. It was then shown that 

{fv{V{^o))-^ol)v' ^0, (121) 

implying V' \\ r°'{V{^o)), contradicting 

Xv{V{^o))V' = 1, (122) 

(which follows from (jll6p ). due to linear degeneracy. 

However, for a linearly degenerate eigenvalue of multiplicity pa > 1, (|117p only guarantees m — 
Lipschitz components, and there are Pa components that we have no regularity information about (indeed at 
a contact with fixed , there is a pa-dimensional submanifold of allowable to which a jump discontinuity 
can occur). Therefore, for the previous argument to work we need differentiability after restriction of a 
function taking values in W" . Saks's theorem does not apply, and functions resembling space filling curves 



show explicit counterexamples. In [1^, the first author proved an analogous result — for any function 



F : D A ^ KP, the set of points x at which 

,\F(x + h) - F{x)\ , , 

lim inf ' ^ , , , ^ = oo (123) 

is of measure zero. This paper also presents an explicit example demonstrating the sharpness of the exponent 
/c/p, and so we cannot hope to obtain differentiability after restriction for an arbitrary function _F : M D 
A — >■ W°' if Pa > 1. Therefore, the different approach of constructing the "intermediate state" in Section 
113.21 is necessary. 



A Regularity of Eigenvalues and Eigenvectors 

The standard implicit function theorem argument for smoothness of simple eigenvalues does not work if there 
are repeated eigenvalues. (We will use the convention that when the eigenvalues are indexed with a subscript 
they are repeated according to their multiplicity, where superscript indices only label the distinct eigenvalues.) 
It is well known (see [28l | ) that the unordered set of eigenvalues (repeated according to multiplicity) of an 
TO X TO matrix is a continuous function of the matrix entries, with the spectrum being an element of C™\ 
where {^aS^^i ^ {/^q}™=i if = l-ia{a) for all a = l..m and some a G S'„. The metric on this quotient 
space is given by 

d({Aa},{/^Q}) = min max \na - K(a)\- (124) 

<tSS„ l<a<m 

However, if the matrices in question are continuous functions of say z E D CZ M'"' such that the eigenvalues 
are real for all z G D, then it is clear we can label 



Ai(z) < X2{z) < ... < A,„(z) (125) 

such that Aq,(z) is a continuous function of z for a = 1, ..,to. 

The smoothness of the eigenvalues and eigenvectors is more delicate when the eigenvalues are not simple 
— in fact there are many counterexamples. However, if the matrices A(z) have the property that each distinct 
X"{z) has constant algebraic multiplicity pa and pa linearly independent eigenvectors for all z G D, then 
around each zq € D there exists a neighborhood D^g 3 zq such that, for a, /3 = 1, .., n, i = 1, ..,pa,j = 1, ..,p^, 

A"(z):L>^R (126) 
r"'*(z) : ^K" (127) 
l""'{z) : D,^ ^R"" (128) 
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are smooth functions satisfying for all z S 

A(z)r"''(z) A"(z)r"'*(z), (129) 

r'\z)A{z) = r'\z)X'{z), (130) 

r^\zy^^{z) = i^Si^,,, (131) 

|r"-*(z)| = 1. (132) 

Moreover, the set of right (and left) eigenvectors is linearly independent for all z, and for each given family 
the right eigenvectors can be taken to be orthonormal. A proof of these statements for a single semisimplc 
eigenvalue of constant multiplicity can be found in [23j . 

Theorem 9 (Nomizu (23] )• ief D C M'^ 6e open and A : D ^ Mm(K) be a smooth mapping such that 
A{z) is diagonalizable for all z. If X is a continuous function on D such that for every z ^ D the value 
X{z) is an eigenvalue of A(z) with the common multiplicity p, then A is smooth. Furthermore, for each zq, 
there exist smooth eigenvectors r^, ...,rP of a neighborhood Dz„ of zq into such that, for each z € D 
r^{z), ...,rP{z) form an orthonormal basis of the eigenspace of A{z) for X{z). 



Zq ; 



Apply this theorem to each eigenvalue, and take the left eigenvectors to be the rows of the inverse matrix 
of the matrix of right eigenvectors to obtain the normalization p3ip . 



B Calculations for Full Euler 

The full Euler equations are given by 



I 



m 

m^p^^ + p 
ranp^^ 



mpp 



) 



n 

mnp~ 
1 



9 



(133) 



y in(m^ + n^)/9 ^ + ne + npp ^ J 



where p is the density, m and n are the horizontal and vertical momentum densities, p is the pressure, 
and e is the specific internal energy. The conserved quantities are p, m, n, and the total energy per 
unit volume pE := j^{nT? + n^) + pe. For convenience, we can calculate the eigenvalues and check for 
genuine nonlinearity/linear degeneracy by considering the matrix pair (/j^,/^) with convenient variables 
W = {p, TO, n, S) where 5 is the specific entropy (the invariance of eigenvalues and nonlinearity properties 
was discussed at the end of Section |9]). It is well known (see Q) that of the five thermodynamic quantities 
p, e, S, p and temperature T, only two are independent, so once two have been chosen the other three can 
be expressed in terms of those two. We have that 



Ppip,S) ==:c2>0, 

ep{p,S) = ^, 
P 

es{p,S) = T, 



(134) 
(135) 
(136) 



where c is the sound speed. The derivatives of the fluxes in terms of the variables W are then 



Iw 



fV 

Jw 



( 







1 







-mnp 

,2 I „2\„-3 



-2 



y —n{m^+n^)p + nc p 





np~ 


mnp~ 



1 

mp~^ 
2np-^ 

2 3„2„-2 1 1^2„-2 







—m'^p^^ + 2mp~^ 
— mnp^^ np~^ 
\ —m{m? + n'^)p^^+mc^p~^ ^iri^ p~'^ + ^n"^ p^^ + e + pp~^ mnp^^ 



U PS 

mp~^ 



>n-p 



m p 



raT + mpsp ^ J 

\ 


VS 

e+pp~^ nT + npsp~^ / 



(137) 



(138) 
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Using Maple, we can compute the eigenvalues and eigenvectors. There are two simple eigenvalues corre- 
sponding to acoustic waves: 



A±(C/) 



mn ± pcyjrn? + - (pc)^ 



— (pc)^ ' 

and a double eigenvalue corresponding to shear waves and entropy jumps: 

m 



(139) 



(140) 



The background state we consider will be of the form U = (1, A/q, 0, E'o), where the units have been 
scaled so that the background density and sound speed are 1, and so the background Mach Number \Mo\ is 
the absolute value of the background horizontal momentum. So that these eigenvalues are real and distinct, 
we must assume |Mo| > 1, which is simply requiring the background state to be supersonic, and it is clear 
that the eigenvalues and eigenvectors are smooth functions on a neighborhood of this U. The eigenvectors 
considered below are also smooth on a neighborhood of U. 



Two linearly independent eigenvectors for A are 



( \ 

m 
n 

V / 



-ps \ 



V 0^ J 



(141) 



As A° is semisimple of constant multiplicity, Boillat's theorem (Theorem [T]) guarantees that A" is linearly 
degenerate, though an easy calculation confirms this as well. For any contact discontinuity corresponding 
to A°, the flow will be tangential to the discontinuity. The r^'^ field corresponds to a shear wave, which is 
a discontinuity in the tangential velocity, and the r"'^ field corresponds to an entropy wave, in which the 
entropy and density are different but in a manner so that the pressure is constant on cither side. 

The acoustic characteristic fields A^ are genuinely nonlinear — to confirm this it suffices to check at the 
background state with n = since genuine nonlinearity is an open condition. From Maple we have that 
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However, the first three entries are precisely the r^ eigenvectors for the isentropic case investigated in [ij], 
A^ are identical in the full and isentropic cases, and the first three W variables we consider in this case are 
precisely the conserved quantities for the isentropic case. Since A± for full Euler are independent of S, it is 
clear that in our context of full Euler A 



w 



is identical to X^r^ from the isentropic case in Section 18 of 



[l3 |. Genuinely nonlinearity was demonstrated for Cp > —1, which is true for commonly used equations of 
state (in particular, any polytropic gas with 7 > 1.) 

As discussed in Chapter II, Section 1.1 of [l3|, taking r]{U) = ~pS, 'ip^{U) ~ —mS, ip'^iU) ~ —nS yields 
an entropy-entropy flux pair. The second law of thermodynamics implies that S* is a strictly concave function 
of and e, which [l7| proves is equivalent to —pS being a strictly convex function of p, m, n, and pE. 
Finally, in Chapter II, Section 2 of 17[, it is shown that the eigenvalues of are 
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(A change of coordinates was introduced that preserves the eigenvalues, since the actual expression for 

2_|_ 2 

is complicated due to the need of differentiating with respect to ™ — h pe.) Therefore, if Mq > 1, eyv is 



positive definite (by Lemma [T]) and the system is hyperbolic in the positive x-dircction and forward sectors 
lie in X > 0. If Afo < —1, eyv is negative definite and the system is hyperbolic in the negative x-direction. 
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